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that an energy influx from bulk Unruh radiation to the dust matter on the brane accelerates 
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I. INTRODUCTION 

An unexpected discovery appeared in 1998, that is, our universe is accelerating rather than 
decelerating This discovery has significant and far-reaching influence on both fundamental 
physics and astronomy. The invisible sector of the universe (or the modified terms to Einstein 
gravity) dominates the present evolution and determines the final destiny of the universe. To 
understand its physical nature, we have to go beyond the standard model. Various explanations 
of this acceleration have been proposed; see for recent reviews with fairly complete lists of 
references of different models. However, the physical nature of the cosmic acceleration remains as 
a most profound problem in sciences. 

In different approaches to the cosmic acceleration problem, braneworld models draw much 
attention in recent years. Inspired by string/M theory the braneworld gravity is introduced for 
several problems in standard model, firstly for the hierarchy problem. The braneworld gravity 
has been extensively studied, for a review, see^. In the braneworld scenario the standard model 
particles are confined to a brane, while gravity can propagate in the whole space time. At the low 
energy region general relativity recovers, while at the high energy region the behavior of gravity is 
strongly modified. This may yield remarkable changes in gravity dynamics with several different 
implications for cosmology, black holes and high energy physics. Among various braneworld models, 
the model proposed by Dvali, Gabadadze and Porrati (DGP) Q] is impressive and a leading one in 
cosmology. In the DGP model, the bulk is a flat Minkowski spacetime, but a reduced gravity term 
appears on a tensionless brane. In this model, gravity appears 4-dimensional at short distances 
but is altered at distance large compared to some freely adjustable crossover scale tq through 
the slow evaporation of the graviton off our 4-dimensional brane world universe into an unseen, 
yet large, 5th dimension. At short distances the 4-dimensional curvature scalar R dominates and 
ensures that gravity looks 4-dimensional. At large distances the 5-dimensional curvature scalar ^^^R 
takes over and gravity spreads into extra dimension. The late-time acceleration is driven by the 
manifestation of the excruciatingly slow leakage of gravity off our four-dimensional world into an 
extra dimension. This offers an alternative explanation for the current acceleration of the universe 
Q]. However, just as LCDM model, it is also suffered from fine-tune problem. In LCDM model, 
A is an unimaginatively tiny constant compared to the expectation of the vacuum if our quantum 
field theory (QFT) is valid till to Planck scale. We know QFT is at least valid at eletro-weak 
scale, whose expectation of zero-point energy is still 10^^ higher than the scale of the cosmological 
constant. To solve the dark energy problem in frame of DGP model, people require Tc ~ Hq^, 



3 



where rc is a constant, and Hq is the present Hubble parameter. Similarly, we can ask why Tc is 
so large that it is approximately equals the Hubble radius. 

Interaction is a universal phenomenon in the physics world. An interaction term between 
bulk and brane has_been invoked as a possible mechanism for the cosmic acceleration in frame of 
braneworld model However, in previous works little attention has been devoted to the physical 
mechanism of such a term. There is a natural physical origin for this interaction which has not been 
adverted before, that is, Unruh effect. In the braneworld scenario a brane is moving in the bulk, 
and generally speaking, its proper acceleration does not vanish. Thus, a brane, qua an observer 
in the bulk space, should perceive the Unruh radiation. The Unruh effect of RS II braneworld is 

■ ■ Q 

mvestigated m [71J. 

Even for a Minkowski bulk, the brane will sense a thermal bath filled with radiations. If the 
temperature is different from the inherent temperature of the brane, energy flow will appear based 
on the most sound principles of thermodynamics. There are several different particles are confined 
to the brane, which are not in thermal equilibrium. Then, what is the inherent temperature of 
the brane in the frame of gravity theory? It will be useful to give a brief account of the previous 
approaches, especially on black hole thermodynamics. 

The relation between gravity theory and thermodynamics is an interesting and profound issue. 
The key quantities bridging the gravity and thermodynamics are temperature and entropy. Tem- 
perature of an ordinary system denotes the the average kinetic energy of microscopic motions of 
a single particle. To gravity the temperature becomes subtle. Since we do not have a complete 
quantum theory of gravity, for general case we can not use the usual way to get the temperature 
of the gravity field. Under this situation people set up some thermodynamics and statistical quan- 
tities of gravity by using semi-quantum (matter field is quantized, but gravity remains classical) 
theory, though the concept of gravitational particle is not clear. The black hole thermodynamics 
(in fact, spacetime thermodynamics, because the physical quantities in black hole thermodynamics 
should be treated as the quantities of the globally asymptomatic flat manifold) is set up in [^i] and 
confirmed by Hawking radiation ^. 



In the initial work [8!], the temperature of the spacetime is suspected by an analogy to the 
ordinary thermodynamics. Shortly, it is recognized that the temperature is just the temperature 
of the radiation emitted from the black hole by a semi-quantum treatment. The method in 
depends on the global structure of the spacetime. The vacuum state is treated as nonvacuum state 
by observers at spacelike infinity, which roots in the fact that vacuum state around the black hole 
is different from the vacuum at spacelike infinity. Along this clue Unruh found that a uniformly 
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accelerating particle detector with proper acceleration A in Minkowski vacuum perceives different 
excited modes exactly in Bose-Einstein distribution with a temperature T = A/2tt [10]. Unruh 
effect helps us to derive the local temperature of a gravitational system. Unruh effect has inherent 
relations to Hawking effect. After a conformal transformation, the thermal particles detected by 
an accelerated detector becomes thermal particles seen by an inertial observer in curved spacetime, 
which is just the Hawking effect in the case of black hole [111 ], 

The laws of black hole mechanics are results of classical Einstein gravity, for qualification of 
laws of thermal dynamics the quantum theory is required. The mathematical form keeps invariant. 
This implies that the classical gravity theory may hide information from quantum theory. This 
possibility is investigated in [2]. The Einstein equation is reproduced from the proportionality of 
entropy and horizon area together with the first law of thermal dynamics, 5Q = TdS, jointing to 
heat, entropy, and temperature, where the temperature is the Unruh temperature to an observer 
just behind a causal Rindler horizon. The entropy is supposed to be proportional to the area of 
this horizon and the heat flow is measured by this observer. The significance of this derivation is 
that the Unruh effect is a result of quantum field theory but the derivation of classical Einstein 
theory depends on it. 

More directly, for a dynamical spacetime, a similar procedure reproduces the Friedmann equa- 
tion. In this case one should apply the first law of thermodynamics to the trapped surface (apparent 
horizon) of an FRW universe and assume the geometric entropy given by a quarter of the apparent 
horizon area and the temperature given by the inverse of the apparent horizon [13]. There are 
several arguments that the apparent horizon should be a causal horizon and is associated with the 
gravitational entropy and Hawking temperature [14]. So, for an expanding universe, the Hawking 
temperature at the apparent horizon should be treated as the inherent temperature of the universe. 
For a brane universe its inherent temperature may be higher or lower than the temperature of the 
Unruh radiation in the bulk, which thus triggers an energy flow between brane and bulk. We will 
study this possibility in frame of DGP braneworld model. 

This paper is organized as follows: In the next section we give the basic construction of DGP 
model and present the Unruh temperature perceived by an inertial observer (inertial respective to 
the brane, which we will explain in detail) on the brane. In section HI, we investigate the condition 
for the temperature to be valid. In section IV, we study geometric temperature of a DGP brane in 
detail and find the thermal equilibrium condition for the Unruh radiation and the brane (evaluated 
by its geometric temperature). In section V, we present a detailed study of the energy exchange 
between bulk Unruh radiation and the brane matter by statistical mechanics. Our conclusion and 
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discussion appear in section VI. 

II. UNRUH EFFECT FOR A DGP BRANE 

In this section we discuss the Unruh effect for a DGP brane. Before studying the Unruh effect 
of a brane, we shah first give a brief review of the gravitational thermodynamics and Unruh effect 
in 4 dimensional theory, especially the reality of this significant effect. 

Unruh effect states that the concept of particle depends on observer. This amazing effect 
generates some puzzles, sometimes even treated as paradoxes. 

First, Unruh's original construction is not consistent because his quantization is not unitarily 
equivalent to the standard construction associated with Minkowski vacuum. Hence some authors 
used mathematically more rigorous methods to solve this problem [l^. 

Second, the temperature of an Unruh particle detector, which is in thermal equilibrium with 

Unruh radiation it experiences in Minkowski vacuum, is higher than the temperature of the vacuum 

to an inertial observer. Hence, does it really emit radiation for an inertial observer, just like an 

accelerated charged particle? It was argued that there was no radiation flux from an Unruh observer 

lla l. Unruh made an almost the same calculation as in tldl, but he found extra terms in the two 

□ 

point correlation function of the field which would contribute to the excitation of a detector IITH . 
It was pointed out that the extra terms were missing in p^, which shape a polarization cloud 
about an Unruh detector |18]. The above analysis and detailed analysis (including non-uniformly 
accelerated observers) in [19] showed that in a 2-dimensional toy model, there is no radiation flux 
from the detector. Recently it was found that there exists a positive radiated power of quantum 



nature emitted by the detector in 4 dimensional space 

It should be noted that the response of an Unruh observer to the Minkowski vacuum is inde- 
pendent of its inner structure: The distribution of the different excited modes perceived by the 
observer depends only on the acceleration of the observer. This implies that the Unruh tempera- 
ture is an inherent property of the quantum field. The detector only plays the role of test particles. 
This effect still exists even without a detector, which can be seen clearly from the derivation by 
Bogolubov transformation. The Unruh effect can be also derived in terms of the spontaneous 
excitation of accelerating atoms 
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Several contemporary research themes, for example, the black hole thermodynamics [24], the 
quantum entanglement state |2^, the Lorent z sy mmetry breaking [2J]. For extensive references of 
the Unruh effect, see a recent review article [27 1. 
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Before studying the Unruh effect for a DGP brane we warm up by a short review the DGP 
brane. The basic construction of DGP model can be written as follows jj], 



S = S'bulk + 'Sbr 



where 



5buik = J d^x^J^g 



and 



'brane 



d'^x^/—g 



-^K^ + LhTa.ne{9al3,tp) 



(1) 



(2) 



(3) 



Here k'^ is the 5-dimensional gravitational constant, ^^^R is the 5-dimensional scalar curvature. 
The induced metric g^i, is defined as 



(4) 



where the lower case of Greeks run from ~ 3, and n^j is the normal to the brane. (/U = 0, 1, 2, 3) 
are the induced 4-dimensional coordinates on the brane, is the trace of extrinsic curvature on 
either side of the brane and ibrane(fi'«/3 5 V') is the effective 4-dimensional Lagrangian, which is given 
by a generic functional of the brane metric ga/3 and matter fields ip on the brane. In this article we 
adopt a mostly negative signature. For DGP model, an induced 4 dimensional Ricci scalar term 
appears in the brane Lagrangian Lbrane, 

,2 



(5) 



where n is the reduced 4 -dimensional Planck mass, R denotes the scalar curvature on the brane, 
and Ljn stands for the Lagrangian of matters on the brane. 

Assuming a Friedmann-Robertson- Walker (FRW) metric on the brane, we can derive the Fried- 
mann equation on the brane (for a more extensive study of the DGP cosmology in which a 



Gauss-Bonnet term a Weyl term appear in the bulk, see 



k 



p 2 2e 



1/2 



(6) 



where p denotes matter energy density on the brane, Vc = k^^u^, denotes the cross radius of DGP 
brane, e = ±1, represents the two branches of DGP model, H is Hubble parameter, k is spatial 
curvature, and a is the scale factor of the brane. The branch e = 1 was treated as an unstable 
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branch based on the hnear order perturbations. However, a recent study shows that ghosts do 
not appear under small fluctuations of an empty background, and conformal sources do not yield 
instabilities either [29 1. 



Now we begin to investigate the Unruh effect for a DGP brane in a Minkowski bulk. We work 
in a bulk-based coordinate system, 

dS^ = df^ - dR^ - R^dnl (7) 

where dS"^ is the 5 dimensional line-element, dQ^ denotes a 3 dimensional sphere. A free falling 
observer on the brane, which is just the comoving observer in the FRW universe, can be described 
in the FRW coordinates, 

X'^ = (t, Xq, Xq, Xg), (8) 

where t = xq is the standard world time function, and three coordinates of the unit 

3-sphere O3. For a comoving observer, they are constants. Hence, clearly its acceleration is zero in 
the brane- viewpoint. However, it accelerates in the cross direction of the brane, or the direction of 
extradimension. The observer moves along a non-geodesic orbit measured by the full 5-dimensional 
metric. This is the exact meaning of "A brane accelerates in the bulk", and the acceleration of the 
brane indicates the acceleration of such comoving observers. 

We assume a free-falling particle on the brane, as a particle detector, moves along the worldline 

x''{t) = {f,R,xl,xlxl), (9) 

where the capital Latin letters run from ~ 4. For a moving brane in the R direction, the velocity 
of a free falling observer on the brane reads 

C/^ = (f,i?, 0,0,0), (10) 

which satisfies the normalization condition, 

UaU^ = 1. (11) 

The induced metric on the brane reads, 

ds'^ = dS^ -U0U = dt^ - a{tfdnl (12) 

where we have labeled the radius coordinate i? by a new symbol a, which is more frequently used 
in cosmology. 
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The unit normal of the brane is vertical to the velocity and spacelike, 

uaU^ = 0, n^n^ = -1. (13) 
The acceleration of the brane is defined as, 

= U^VbU^- (14) 

We consider the component of the extrinsic curvature along the velocity, 

Ktt = KabU^U^ = -ncA^. (15) 

By using the junction condition across the brane 

[K^,u] = -K^{sf,v - '^g^,vs), (16) 

where [K^^] denotes the difference of the extrinsic curvatures of the two sides of the brane, s^j/ 
represents the energy-momentum confined to the brane, we derive the amplitude of the acceleration 
of brane js^ , 

A = |(5)5CD^^^^|'^' = y \2pe + 3pe| . (17) 
Here pe and pe are the effective energy density and pressure of the brane, which are defined as 

Pe = sl (18) 
Pe = -4, (19) 

where i is an arbitrary spatial index and we need not sum over the index i. And s is the effective 
energy momentum confined to the brane, 

where -Lbrane is given by Due to the above equation, the energy momentum s^i, includes the 
contribution of the matter term L^, as well as the induced curvature term R on the brane. That 
is the reason why we call the density and pressure in (llSp and ()19p effective density and pressure. 
We stress that the effective density pe is different from the local density p of the matter confined 
to the brane in ([6]), which does not include the geometric effect of the induced Ricci term R. 

We see that, generally speaking, A does not vanish. So the brane may perceive Unruh-type 
radiation in the bulk. To discuss Unruh radiations seen by the observer on the brane, we introduce 
a 5-dimensional scalar field in the bulk, which is in its Minkowski vacuum state, 

c^ = ldAcl>dU-W^^, (21) 
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where m denotes (p^s mass parameter. Since the scalar is in its vacuum state, classically it does no 
work on the brane and bulk dynamics. 

Several different methods have been proposed to derive Unruh effect since Unruh's original 
work [11]. Here we use the Green function method. We consider a detector minimally coupled to 



a 5-dimensional scalar field 6 in the bulk 



The Lagrangian cm{t)(f)[x{t)] describes the interaction between detector and field, where c is 



the coupling constant and m denotes the moment operator of the detector. It is shown in [111 ] that, 
for a small number c, the probability amplitude of the transition from the ground state \Eq > of 
a particle detector coupled to a scalar field in its vacuum state ]0a/ > to an excited state \E >, 
where M stands for Minkowski, reads 

c^Y,\<E\m\OM,Eo>\^HE-Eo), (22) 

E 

where the detector response function J-{E) is defined as 

/oo poo 
dt / dt'G^{x{t),x{t'))e~'^^-'^°'^^'-''\ (23) 
-oo J —oo 

and the Wightman Green function G^{x,x') is defined as 

G+ =< Om\HxMx')\Om > . (24) 

We note that the term ] < E\m\OM, Eq > ]^ depends on the structure of the particle detector, but 
J^{E — Eq) does not, which reflects the inherent properties of quantum fields. It is clear that for a 
general trajectory the response function is not zero. 

For a uniformly accelerated trajectory the Wightman Green function becomes a function of 
At = t — t' . Since the detector will detect infinite particles in its whole history, the response 
function (j23p is not well-defined. Under this situation it will make sense to consider the unit 
response function, which describes in unit time interval, 

/oo 
dAte-^(-^-^")^*G+(At). (25) 
-OO 

constant, we just set Xq = Xq = Xq = after a coordinates transformation. For 
a uniformly accelerated observer in the cross-brane direction, the coordinates in ([9]) read, 

190 - sinh(tyl) ^ cosh(tA) , , 

x\ = xl = xl = ^, T= ^ \ R= (26) 



where A is a constant, denoting the magnitude of acceleration of the detector. 



(p^ directly with proper boundary condition (corresponding to proper contour [ll|] ) , we obtain its 



len, integrating 
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concrete form in McDonald function (Bessel function with imaginary arguments), 

(At) = e^^"/"(m^)^/^ ,,,,, i^3/2 ( ^ sinh( AAt/2 - ie)) , (27) 

^ ^ 167r5/2(sinh(^At/2-ie))^/2 -^'^ \ A ^ ' ^ ' 

where, as usual, e is a small positive number. The Wightman Green function for massless mode 

reads, 

D+(At) = lim G+(At) = — r- — rmr-, (28) 

^ ^ m^o ^ ^ 647r2 (sinh(^At/2-ie))3 ^ 

Generally speaking, Unruh effect is a very weak effect. For example, in 4 dimensional Minkowski 

space, it needs about an acceleration 10^^ meter per second^ to increase temperature 1 K, hence 

various massive modes are difficult to be excited. The mode with zero mass is the easiest mode to 

be excited. Therefore, we consider excitations of massless modes. Substituting D^{At) into (j25p . 

one obtains the unit response function with a contour closed at lower half plane of complex At by 

using the Jordan's lemma, 

- - 2,2^ e^AE-E,)IA + 1 • (29) 
This is the response function for a particle detector confined to the FRW brane, and comoving to 
the FRW brane, which describes photon gas system at temperature Aj^ix. We see that a Fermi- 
Dirac factor appears, which implies that the particles of Rindler radiation behaves as Fermions, 
though all the excited modes we integrated to derive the response function are Bosonic. It is not 
a phenomenon completely new. In 1986, Unruh pointed out that the Fermi-Dirac factor would 
appear in the response function for an accelerated monopole of a massless field in an odd number 
of space dimension, arising from integration over all modes for a scalar field 1^. 

In the above text of this section, only a brane with positive spatial curvature is discussed. A 
spatially flat brane can be treated as a limit when K is large enough. To the case of a negative 
spatial curvature, we need to replace the 3-sphere by a 3- hyperbola in ([7]). A 5-dimensional 
Minkowski space sliced with space-like 3-hyperbola can be written as, 

dS^ = -df2 + dR^ - R^dBl. (30) 

Then the following discussions exactly follow the case of positive spatial curvature. The resulting 
distribution function (1291) is still valid. 



III. QUASI-STATIONARY ACCELERATION STAGE 



Up to now our analysis is limited to uniformly accelerated detectors. But from (jl7p we see that 
the acceleration of the brane is not a constant. Physically, the formula ()29p remains valid when 
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the acceleration varies slowly. Mathematically we can estimate the time over which the constant 
acceleration approximation is valid by imposing that the variations in the acceleration are small, 
e.g., expanding the acceleration around some time Iq to the first order t — to: 

dA 

A{t) = A{to) + {t-to)—\t=to (31) 

Here, we adopt the cosmic time (proper time) t rather than the proper time coordinate time T of the 
detector since the cosmic time is directly related to our measurement of cosmological parameters. 
For evaluating the variation of the acceleration, define characteristic time 

dA/dt\t=t, ' ^^'^ 
which is a function of Iq, i.e., it varies with the evolution of the universe. Our condition for 
"acceleration varies slowly" requires that the time scale of a physical process we concerned is much 
less than the characteristic time, 

ti <C tc. (33) 

The time scale of a physical process ti must be shorter than the age of the universe. So if the 
characteristic time tc is much larger than the time scale of the universe, we can safely treat the 
acceleration as a constant for any physical processes. 

Here we present some examples to explain this condition. We work in a frame of "rcCDM" 
model, that is, a DGP brane universe filled with dust, whose density function p in Friedmann 
equation only includes pressureless matter, with about 30% of the critical density. 

First, we consider the inflationary phase at the early universe. At such a high energy scale, the 
infrared correction of DGP model to the standard general relativity, which becomes important at 
late time universe, can be omitted safely. So the effective energy density and pressure in (jl7p are 
just the ordinary density and pressure of the universe. At the whole inflationary phase, the density 
of the universe is approximately constant, and the pressure p = —p. Hence, 

^ = --3^-0, (34) 
dt 6 dt ' ^ ^ 

that is, from ()32p . the characteristic time is very long. In this case, the Unruh temperature is well 

defined. 

Second, the universe enters a radiation dominated phase after the inflation. At this stage, the 
DGP-correction to the standard model is still tiny. Thus, we adopt the the same approximation 
as above, pe = p, Pe = P- For a radiation dominated universe, we have 

P = ^P- (35) 
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The continuity equation reads, 

dp 



^^+3H{p + p) = 0. (36) 
Substitute ([35ll and ((36]) into (f32]) . we reach, 

from which we see that the characteristic time becomes longer when the energy density becomes 
lower. For instance, when p = (O.lMev)^, we get tc ~ Is, which is much shorter than the age of 
the universe at that time. 

And then the universe is diluted to be thinner and thinner. The DGP correction becomes more 
and more important. Under this situation we have to introduce the correction terms. By using 
([201) (see also [31]), we derive 

Pe = p-p\3X), (38) 
p,=p + p\^+X), (39) 

where X is defined as 

X^H^ + k/a^. (40) 
Substituting to the Friedmann equation Q, we obtain 

As for standard cosmology, ^ can be obtained from the Friedmann equation and the continuity 
equation . Here, similarly, from Friedmann equation ([6]) and the continuity equation (|75p . we derive 



'^ = X-±,^±P^. (42) 



Then substituting ([38]) , ([39|) and ([T7|) into ([32]) , we arrive at the following rather complicated 
from, 



1 

tc= - 

3 



3p^X ^. e 



2 + ^ 3(1 + 



k 



(X--2)-VV-\ (43) 



where 
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Substituting the current parameter of the universe p = (1.8 x 10~^ev)^ and we suppose that the 
universe is spatially flat, then we obtain, either for the negative or the positive branch, tc ~ 10^''s, 
which is at the same scale of the age of the universe. So the constant acceleration is a perfect 
approximation under this condition. 

Then, at least of late time universe, we can safely use (I29p . which is a thermal distribution at 
a temperature of 

^ = (^^' 

Associating with pT|) . we get 

:|2/9e + 3pe|. (46) 



127r' 



IV. COMPARED TO GEOMETRIC TEMPERATURE OF THE BRANE 

We see that the temperature of Unruh radiation the brane perceives, as displayed in (j46p . 
only relates to the energy density and pressure of the brane for a given DGP model. Then, a 
natural question emerges: Is the Unruh radiation hotter, colder or at the same temperature to 
the brane? If there is energy exchange between bulk and brane, we also need the temperature of 
the brane to decide the direction of the energy flux. But there are several different particles on 
the brane. They had gone out of thermal equilibrium long before. Hence we'd better to find the 
characteristic temperature of the braneworld independent of its detailed microscopic construction. 
As we have pointed out Unruh radiation is a geometric feature of the space, which is unrelated to 
the construction of the detector. Just as well we need the characteristic temperature of the brane. 

We know that the formulae of black hole entropy and temperature have a certain universality in 
the sense that the horizon area and surface gravity are purely geometric quantities determined by 
the space geometry, once Einstein equation determines the space geometry. As we have mentioned 
above, Einstein equation can be reproduced by thermal dynamics considerations. As for the case of 
cosmology, applying the first law of thermodynamics to the apparent horizon of an FRW universe 
and assuming the geometric entropy given by a quarter of the apparent horizon area and the 
temperature given by the inverse of the apparent horizon, the Friedmann equation can be derived. 
This celebrated result implies that the inverse number of the apparent horizon is the geometric 
temperature of the universe, which is independent of the microscopic structures of the particles 
confined to the brane. The apparent horizon in the dynamical universe is a marginally trapped 
surface with vanishing expansion. Straightforward calculation yields the radius of the apparent 
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horizon 



(47) 



where X is defined in (|40p . And then the geometric temperature of the brane reads 

-1 



1 



2-K 27r 2tt V 3^*2 r2 \?,IJ? 
Substituting (glD into ^ and ([MI) and ([Ml) into (gSD, we arrive at 



1/2 



(48) 



(49) 



127r^2 

where X is given by (j4ip . Therefore we obtain the exphcit form temperature of the brane, which 
is determined by its energy density and pressure. By contrast the geometric temperature only 
depends on the energy density, which has no relation to the pressure. We see that both Unruh 
temperature (|19|) and geometric temperature ([l8j) are unrelated to the spatial curvature of the 
brane. 

To investigate further the evolution of the Unruh temperature and the geometric temperature 
during the history of the universe, we write them in dimensionless form, 



Ho ~ 2^Ho 



1 

2^ 



oc^^Yfi I 2^^j-'^ I 2c -^/^^-j^^ i^^x^^ffL I ^^'ff^ ^ '''' 



1/2 



(50) 



T 



A 



[l + ^w/2)+ ^ 



27r./n 



"^^"^ :(l + 3w;/2)+ ^ 



—2 - 3J^ma;(l + u;) — 



27r^/J^ 



47rA/0^ 



47r2r2 



3r2ma;(l + ui) 



1 



;;,i) 



where Hq is the present value of the Hubble parameter, w is the parameter of state equation of 
the matter confined to the brane, x = p/po, and Qm, are defined as 

Po 



3p^Hi ' 
1 



(52) 
(53) 



Here po denotes the present density. First, we study various limits of the two types of the tem- 
peratures. When Tc ^ oo, ^Irc ~^ 0, we expect that DGP brane theory reduces to standard 
4-dimensional cosmology. It is right in the case for geometric temperature of the brane, 

VP 



lim r' - ^ . 



(54) 
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But for Unruh temperature, it is a different situation. Though all dynamical effect of the 5th 
dimension vanishes when Tc — > oo, the Unruh temperature does not vanish. One can check 



1 ^ \fP 

lim T = — |5 + 3t(j| — . 
'"c^oo 2 27rV3^ 



(55) 



This result implies that even when the gravitational effect of the extra dimension vanishes, the 
quantum effect is saved. The other important limit is the limit when the matter on the brane is 
infinitely diluted, i.e., p — > 0, 



One may conclude 



lim T 



27r Tc 
1 1 + e 



lim T 

p^o 27r r, 



(56) 
(57) 



for either e = 1 or e = — 1. However 



for e = 1, and 



lim T' = lim T, 
p^O p^O 



T' 

lim — = 1, 

P^o T 



, T' 1 
lim — = ■; r, 

P^o T 2 + 3t(; ' 



(58) 



(59) 



(60) 



for e = — 1, since the "speeds" are different when T and T' go to zero. 

An important case for cosmology is a universe filled with dust matter, i.e., w = 0. In this case, 
ED and (09]) become. 



T' = Z + eY, 



(61) 



and 



where we set 



22Y + eZ 



(62) 



and 



ZTT 



(63) 



(64) 
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For a reasonable cosmological model, we have y > 0, Z > and Z > Y. The calculation to take 
absolute value is difficult to deal with. Hence we consider different cases in which we can determine 
the sign in the absolute calculation. 
First, in the branch e = 1, T becomes 

because every term is larger than zero. Then 

Therefore, the Unruh radiation perceived by a brane is warmer than the brane measured by the 
geometric temperature. 

Second, we consider the branch e = —1. The Unruh temperature in this case is much more 
complicated. First, at the point 2Y = Z, or equally 30,rc = xQm, there is a singularity T ^ oo 
in (j62p . which means the particle numbers are equal at every energy level. This is valid only for 
a system endowed with finite energy levels, which indicates our theory must be cut off at some 
energy level. On the right hand of this singularity 2Y — Z > 0, T can be decomposed in the form 

{Y - Z){Y - 5Z) 



2{2Y - Z) 



(67) 



where Y — Z < Q,Y — bZ < and 2Y — Z > 0, hence 2^Y^Z)^^ ^ ^- remove the absolute 

sign directly and compare T and T' by 

(y-z)(5y-7Z) 

2{2Y-Z) ^''^ 



since Y - Z < 0, 5Y - 7Z < and 2Y - Z > 0. On the left hand of this singularity 2Y - Z <0, 

(Y 

similar to the case 2Y — Z > 0, we can prove — 
by inserting a minus sign and, T — T' becomes 



similar to the case 2Y — Z > 0, we can prove — 2(2Y-Z) — ^ ^' Hence we remove the absolute sign 



72 _ T^2 

Thus we prove that the Unruh radiation perceived by a dust dominated brane is warmer than the 
brane measured by the geometric temperature all the time. 

The other point deserving to be noted is that T' will be always equal to T if the parameter of 
state equation of the matter is confined to brane w = —1. In this case. 
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FIG. 2: u = 2ttT/Ho as a function of u = 2'kT' /Hq. In this figure we also consider quintessence-like, in 
which w = —1/2 and we set il„i — 0.3. (a) The branch e = 1. (b) The branch e = —1. 

For general case the expressions of T and T' are rather complicated. Hence we plot two figures 
to give their visual profiles. Fig. 1 illustrates T'/T as a function of x and 17^^. As we have pointed 
out, when /) — > 0, T'/T — > 1 for the branch e = 1; while when p ^ 0, T'/T 1/|2 + 3w\ = 1/2 
for the branch e = 1. Fig. 2 directly displays T as a function of T' , in which the interval of the 
argument x is x G (0,3). Recalling x = p/po, we see that we consider the temperatures in some 
low redshift region and in the future in figure 2. It is clear that T ^ when T' ^ 0, as shown by 
(I56p and ()57p for the branch e = 1. Our numerical result also shows an interesting property of the 
branch e = — 1: T is almost a linear function of T' . 
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V. UNRUH RADIATION ACCELERATES THE UNIVERSE 

With the detailed studies in the last sections, we see that the Unruh temperature does not 
vanish even at the limit Tc — > oo, which is the condition for the vanishing of the gravitational 
effect of the 5th dimension. Furthermore, the Unruh temperature is always higher than geometric 
temperature for a dust dominated brane in the whole history of the universe for both of the two 
branches and for all three types of spatial curvature. So an influx from bulk to brane will appear 
if brane interacts with bulk based on the most sound principle of thermodynamics. 

Although this difference in temperature indicates an energy flux, it offers no hint of the form 
of the interaction term. Under this situation, we explore the microscopic mechanism based on 
statistical mechanics and particle physics to derive the interaction term between Unruh radiation 
and dark matter confined to the brane. 

In section II, we have considered a detector coupling to a scalar field = ^dA(l>d^<f> — ^m'^(jP, 
which satisfies 5 dimensional Klein-Gordon equation. At the massless limit, it obeys 

□(^V = 0. (71) 
By using Green function method, we derive the temperature of the Unruh radiation, 

T=A, (72) 

where A is given by pzp . 

To investigate the microscopic mechanism of interaction between brane and bulk, one should 
decompose the 5-dimensional modes, which satisfies (|7ip . into modes along the brane and modes 



transverse to the brane. Following [3^, we call them type I modes and type II modes, respectively. 
In a pseudo-Euclidean coordinates system, type I mode reads, (j)^-^^ = e"*'^''^''', where /c^ denotes 
the momentum of the mode in different directions, and /i = 0,1,2,3. Type II modes reads, 
(^(^^) = ^-^k^x'^-lkyy ^ where y is the coordinate of the extra dimension, and we require ky ^ 0. An 



^ that Type II 



arbitrary 5 dimensional mode must be either type I or type II. It was shown in 
modes of the bulk fluctuations do not interact with the brane, for which the brane is effectively a 
mirror. Hence only type I mode can exchange energy with matter confined to the brane. We show 
a schematic plan of type I and type II modes in Fig. 3. 

By a standard process of canonical quantization, type I modes become photons on the brane. 
Here photons mean quanta without mass, which do not necessarily satisfy Maxwell equations. In 
fact we do not even need to study its commutativity in detail. Whether they obey Fermi-Dirac 
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Modes transverse to the brane 



FIG. 3: A sketch of type I and type II modes. Only type I mode can exchange energy with particles confined 
to the brane, while the brane behaves as a perfect mirror for type II mode (for a color version on line). 

statistics or Bose-Einstein statistics the number density of the photon gas is proportional to the 
cubic of its temperature, 



where is just the temperature of the Unruh radiation (j45p . The particles confined to the brane 
are immerged in the thermal bath of these photons. Therefore, based on statistical mechanics 
the reaction rate T between dark matter particles and the type I photons is proportional to the 
number densities of photons and dark matter particles, the relative velocity of the photons and 
matter particles, and the scattering cross section F oc n^^^n^avprn- Here, denotes the number 
density of dark matter particle, stands for the relative velocity, which is a constant, and a 
represents the scattering cross section. In a low energy region, the internal freedoms can not be 
excited. This paper concentrates on the late time universe, hence the cross section is effectively 
constant. Therefore, the reaction rate can be written as an equation by inserting a constant b, 



where we have used ([73]) . (|i5]) and p. Here, we assume there is only pressureless dark matter 

on the brane, which interacts with the Unruh photons. Therefore, the continuity equation of the 
brane becomes, 




(73) 



r = bpA^, 



(74) 



p + 3Hp = T. 



(75) 
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As we explained before, the DGP braneworld model also suffers from fine-tune problem, which 
says why the three RHS terms of ([6]) are at the same scale. To overcome this hurdle, we consider 
a limiting DGP model, that is, Tc ^ H^^. In such a model, the fine-tuned problem is evaded, 
and at the same time the gravitational effect of the 5th dimension is no longer responsible for the 
present acceleration. We have seen that the Unruh effect does not vanish even when Tc ^ oo. We 
will prove that under the situation Tc ^ ^^o"^ the interaction between dark matter on the brane 
and Unruh radiation can drive the observed acceleration of the universe. 

In a limiting DGP model , substituting the Friedmann equation ^ into (j75p . we derive 

p k 



-1/2 



(76) 



a p 
a 6/x2 

The acceleration A of the brane in the bulk does not directly depend on energy density of the 
matter on the brane p, but through the effective density pe and and effective pressure pe in (jl7p . 
which are presented in ([38]) and ([39|) . respectively. 
Associating (HZ]), dMI), ([39]), and ([76]) we derive 

The formula (j77p is general for all the three cases of curvatures. In the case of a spatially flat 
brane, it takes a simple form, 

A-^ = C^p, (78) 

where C = K^pb. Interestingly, we see that the acceleration of the brane A inversely correlates to 
the energy density of the brane, which is completely different from the first sight at equation (jl7p . 
Therefore, we expect the effect of the bulk Unruh radiation is negligible in the early time for a dust 
dominated limiting DGP braneworld. Only in some low energy region the bulk Unruh radiation 
becomes important. 

Before presenting the exact cosmic solution, we study some qualitative side of the Firedmann 
equation ([6]) and continuity equation ([75]) to see whether our model is stable. We consider a 
spatially flat universe, in which the stagnation point dwells at 5H = 0, or equivalently 5p = 0. 
From the continuity equation ([75]) . we derive the energy density ps at the stagnation point, 

p5/4 = ft^c-^/V^/s. (79) 

To investigate the stability of the cosmic fluid in the neighbourhood of the stagnation point, we 
impose a perturbation to the continuity equation, 

(80) 



m- = 6p (^Ibp-y^c-y^ - j 
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At the stagnation point, p = ps, hence 

m-\,=,^ = -^-^pl'Hp, (81) 

which means it is a stable point. Now we consider the deceleration parameter, which is one the 
most significant parameters from the viewpoint of observations. Here the deceleration parameter 
reads, 

where A is given by (j77p . In the case of a spatially flat universe, it degenerates to 

q = \(l- VSpbC-'/'p-'/') . (83) 

This equation clearly shows that in the early universe q — > 1/2, hence the universe behaves 
as dust dominated one, and with the decreasing of energy density the deceleration param- 
eter becomes smaller. Finally at the stagnation point the deceleration parameter ceases at 
q = ^ (^1 — \/3pbC~^^^ps ^^^^ = —1, which implies that the universe enters a de Sitter phase. 

Though, for an arbitrary spatial curvature the analytical solution does not exist, we flnd an 
exact solution for a spatially flat universe driven by bulk Unruh radiation. For a spatially flat 
universe, associating Priedmann equation ([6]) with the continuity equation ()75p we obtain. 



+ ^ ) , (84) 



where ci is an integration constant, and L is deflned by L = y/?>phC~^/'^ . When ci > 0, from 
this exact solution, the universe behaves as dust dominated one in a high energy region (a small 
enough), and becomes de Sitter universe in a low energy region (a large enough), which is exactly 
the same as we concluded before from behaviors of the deceleration parameter in the history of the 
universe. When ci < 0, there is a bounce in the early universe, and the universe also enters a de 
Sitter phase in a low energy region. Since our investigations for the acceleration driven by Unruh 
radiation concentrate on the low energy region, it is only a toy model in the early universe. 

With the exact solution (j84p and dG]), we obtain an analytical expression of cosmic time t as 
function of the scale factor a by using a hypergeometric function. 



V^p 32/5 '5' 5' 5' 3cr 
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where C2 is an integration constant, and F denotes Gauss hypergeometric function. The require- 
ment 

hm t = 0, (86) 

a— >0 

yields C2 = 0. 

To visualize the physical meaning of this solution, we demand the various limits of it. In the 
high energy limit (small a), expanding (ISSp in series around a = 0, 



a'/^ 



21/4 + 0(„17/2)^ (87) 



3cT 315cl/^ 

where if we only keep the dominated term we just obtain a ~ as we expected. In a low energy 
region when a is large enough, the series becomes, 

, L , 15 ,2 ■ 

-7 + ln(-) + -lna-^(-) 



VSfi 5 33/5L2/5 

where 7 is the Euler constant, and tp is the digamma function. If we only take the leading term, 
we obtain a = Me^*, where M = exp{^[7 - ln(L/3c) + V(2/5)]}, N = L^/^S^/^iSV^fiy^, which 
is just a de Sitter space. 

Thus, we complete a cosmic solution of the limiting DGP braneworld model, where the universe 
is self-accelerated through the bulk Unruh radiation perceived by the brane. From the viewpoint of 
observations our model belongs to the class of unified dark energy model, that is, there is only one 
component (usually dust) in the universe, but due to different reasons it evolves in a non-standard 
way. Therefore, it can drive the present acceleration of the universe. Such models also have been 
phenomenologically investigated in, for example, js^. 

Finally, we would like present some preliminary discussions to confront the observations. First, 
the parameter b is critical to our model, which encloses all the undetermined information of inter- 
action between dark matter and Unruh radiation. Here we present a preliminary estimation of its 
value from the deceleration parameter (|83|) in a spatially flat universe. From various observations 
the present value of deceleration parameter q ~ —0.5, therefore. 



ircHor'^' i -4-, ) ~ 1, (89) 



-5/4 



V2 

where ^^2jji = 1 in our model. We see that a larger cross radius needs a smaller coupling constant 
b. And we have set VcHq » 1 in previous constructions. Hence the dimensionless coupling constant 
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bHQ is a tiny number, which eludes our laboratory experiments even if the dark matter particles 
were found. Alternatively, in future work we expect astronomical observations to fit our model and 
hence to determine the value of b, which is also helpful to constrain the cross radius in DGP 
model. 

An exotic matter with negative pressure, call dark energy, is frequently introduced to explain 
the cosmic acceleration in frame of general relativity. To explain observed accelerated expansion, 
we calculate the equation of state w of the effective "dark energy" caused by the induced Ricci term 
and energy influx from the bulk Unruh radiation by comparing the modified Priedmann equation 
in the brane world scenario and the standard Friedmann equation in general relativity, since almost 
all observed properties of dark energy are "derived" in general relativity. The Friedmann equation 
in the four dimensional general relativity can be written as 



H^ + ^ = :X{P + Pde), (90) 

CI oLt 



where the first term of RHS of the above equation represents the dust matter and the second term 
stands for the dark energy. Generally speaking the Bianchi identity requires, 

^+3ii'(pde+Pde) =0, (91) 

we can then express the equation of state for the dark energy as 

Pde . Idlnpde ,„„s 

Wde = — = -1 - -— . (92) 

Pde o dina 



2 2e ( p 1 ^ 



Comparing (f90l) and ([6]), we derive 

Note that in our model there is no exotic matter confined to the brane. pde is in fact geometric 
and interaction effect. We call pde equivalent or virtual density of dark energy. Various evidences, 
which are independent to cosmological models, implies the existence of dark matter with present 
density about $7^ = 0.2 ~ 0.4 



34l |. Maybe more or less, but the density of dark matter does reach 
the density to flat the space. In our model, the geometric contribution of can be very small. So 
generally speaking, we need a curvature term. Substituting ([93]) into ([92|) . and recalling that the 



matter confined to the brane is pressureless, we obtain. 



/ 1 \ -1/2 m / 1 \ r, n -1 

e p p ' 



Wde = -1 + - „ n,, TT^ + 



3 3p'^H ySp"^ rl 



(94) 
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where p is given by ([75]) . and T can be calculated by ([7^ . Far a large Tc, after complicated but 
straightforward algebraic calculating, we deduce, 

Wde = -\ + Qa-3/2(l + (17^(1 + z)3 + 17,(1 + zff'^ , (95) 



where Q is defined as. 



Q = ]^(r,H„)-^l\bHl)-^l\ (96) 

and z is the redshift. It is clear that when Qilm^"^ {^m + ^kY^^ > 0.5, the virtual dark energy 
behaves as quintessence; while when QQ.m^'^{Q.m + ^fe)^^^ < 0.5, it behaves as phantom. We 
stress for the second time that the dark energy in this model is only some virtual, not actual stuff. 
The most sensible quantity in observation is the deceleration parameter which is given in (j82p . 
Substituting (i77|) into (182]) . we reach. 



where we also supposed a large re- Note that for both q and tf^ei e vanishes at the large limit. 
This reasonable, since the two branches have no difference at such a limi.t 

As numerical examples, we plot a figure for q for two sets of parameters in fig HI This figure 
illuminates that the universe is accelerating and q goes to 0.5 quickly at high redshift region, ie, 
the universe becomes a dust one. 



VI. CONCLUSION 



Our scenario is that the brane is accelerating in the bulk, such that it is bathed in the Unruh 
radiation in the bulk. The temperature of the Unruh radiation is higher than the brane, which 
yields energy flux from the Unruh radiation to the brane matter, which accelerates the brane 
universe. 

In this article we study the cosmology case of a DGP braneworld model, that is, an FRW 
brane imbedded in the bulk spacetime. We first investigate the Unruh temperature that the brane 
perceives and its "own" temperature. Generally speaking the brane accelerates in the bulk. Hence 
the particles confined to the brane and inertial on the brane should also perceive Unruh radiation 
in the bulk. In this paper we show that for a DGP brane in a Minkowski bulk it is just the case. 
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1.5 



2.0 



FIG. 4: The evolution of deceleration parameter with respect to z. For both of the two curves, = 0.3, 
rife = 0.7. Q — 0.4 on the solid curve, and Q = 0.3 on the dashed curve. 



We investigate the case of the brane, as a particle detector, coupled to a massless scalar field in 
the bulk. As for a point detector, the temperature of the radiation perceived is proportional to the 
acceleration of the brane. 

As for the brane's own temperature, we should seek a characteristic temperature of the brane. 
The temperature can reveal the inherent property of gravity. The Priedmann equation can be 
reproduced by applying the first law of thermal dynamics to the apparent horizon. Thus we take 
the temperature of the apparent horizon as the characteristic temperature of the brane, which we 
called geometric temperature. 

We find that generally speaking the Unruh temperature the brane perceived and geometric 
temperature of the brane are not equal. We compare these two temperatures in various cases. 
We find that for a dust dominated brane, the temperature of the Unruh radiation perceived by a 
brane is always higher than the geometric temperature of the brane, either in the branch e = 1 or 
e = —1, no matter what the value of the cross radius Tc and the spatial curvature of the brane. 
So, generally speaking the brane and the Unruh radiation it perceived can not reach to thermal 
equilibrium. 

As we pointed out before, the validity of Unruh radiation has been confirmed in detail 
Therefore if energy exchange is allowed between bulk and matter brane, an energy flux between 
the bulk radiation and the matter confined to the brane will come forth, which may accelerate our 
universe. In all cases p > 0, the Unruh temperature of the bulk radiation T is higher than than 
the geometric temperature T' for a dust dominated brane, which means an energy influx to the 
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brane can appear. So we study this possibility in section V. 

The DGP braneworld model seems to be a hopeful candidate to explain the cosmological ac- 
celeration. But, as we pointed out before, it also suffers from fine-tune problem. To evade this 
problem, we consider the limiting DGP model, that is, rc > Hq^, which means the gravitational 
effect of the 5th dimension is negligible. Under this condition, the cosmological acceleration does 
not happen in original DGP model under this situation if there is only dust on the brane. By 
considering the possible energy influx from the bulk Unruh radiation to the brane induced by 
the temperature grads, we find the universe can accelerate through interacting the bulk Unruh 
radiation in a limiting DGP model, even the brane is dust dominated. Differently from previous 
works on the brane-bulk interaction, we find the interaction form through careful studies on the 
microscopic mechanism of interaction between brane and bulk. It is shown that the interaction 
term can be settled up to a constant factor b. Based on these constructions we find the acceleration 
of the brane A is inversely correlated with the energy density of the brane for a dust dominated 
limiting DGP brane. Finally we derive an exact solution for a spatially flat model. This solution 
shows clearly that the universe behaves as a dust dominated one at early time and enters a de 
Sitter phase at late time, which is consistent with observations. We also show the de Sitter phase 
is stable. 

In this paper only massless Unruh mode is considered. Although it is the most important mode 
in low energy region, the massive mode also deserves to study further for a full Unruh effect in a 
high energy region. Also, in the enough high energy region, the universe is radiation dominated, 
and at the same time the scattering cross section between dark matter particles and Unruh photons 
becomes temperature-dependent, since the internal freedoms can be excited. Under this situation 
the Unruh effect may be important again, which deserves to investigate further. 

At the observation side, we should not only take a special set of parameters to show the property 
of this model, but constrain the parameters, especially b by various observations in the future. 
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[35] One often adds a conformal coupling term to gravitational field in the action Leon = — But 
there are some subtles in this couple, because there exists two Ricci scalars in the DGP model and 
two types of conformal transformations. Wether the Lagrangian is conformally invariant depends on 
definition. Here we just consider the minimal coupling case. 



